N95-13949 


DECOHERENCE AND DISSIPATION FOR A QUANTUM 
SYSTEM COUPLED TO A LOCAL ENVIRONMENT 

Michael R. Gallis 

Department of Physics, Penn State University/Schuylkill Campus 
Schuylkill Haven, Pennsylvania 17972 
Intemet:mrg3@psuvm.psu.edu 


Abstract 

Decoherence and dissipation in quantum systems has been studied extensively in the 
context of Quantum Brownian Motion. Effective decoherence in coarse grained quantum 
systems has been a central issue in recent efforts by Zurek and by Hartle and Gell-Mann 
to address the Quantum Measurement Problem. Although these models can yield very 
general classical phenomenology, they are incapable of reproducing relevant characteristics 
expected of a local environment on a quantum system, such as the characteristic dependence 
of decoherence on environment spatial correlations. I discuss the characteristics of Quantum 
Brownian Motion in a local environment by examining aspects of first principle calculations 
and by the construction of phenomenological models. Effective quantum Langevin equations 
and master equations are presented in a variety of representations. Comparisons are made 
with standard results such as the Caldeira-Leggett master equation. 


1 Introduction and Motivation 


Decoherence via coarse graining has been studied in the context of quantum measurement theory 
by Zurek[1] and by Hartle and Gell-Mann[2] as a mechanism which leads to the emergence of clas¬ 
sical properties. Recent efforts have focused on the decoherence effects of a heat bath, which has 
also been examined in detail in the study of quantum brownian motion. Decoherence is identified 
as the (effective) suppression of interference terms in the density operator ( p(x , x '), x ^ x'). It has 
been pointed out that most of the models which have been considered are somewhat simplistic and 
cannot reproduce the phenomenological features expected of a system which interacts locally with 
a homogeneous and isotropic environment [3], In this paper I describe the perceived shortcomings 
of existing models and illustrate the construction of a phenomenological quantum master equation 
which contains many features expected from local coupling to a homogeneous environment[4]. 

Although decoherence is the most interesting feature of the effects of a heat bath, dissipation 
(and other effects) also generally appear in the dynamics of the density operator: 


dp(x, x'\ t) 

di 


= Hamiltonian terms -I- Dissipation terms + • • • — g(x,x')p(x, x';t). 


(1) 


The decoherence term appears as a (spatially dependent) decay term in the evolution equation, 
and can be understood in terms of effective fluctuating forces, or potentials: [5, 6] 

9(x,y) = (^ 2 )(c(x-,x) + c(y; y) - 2 c(x;y), (V(x,t)V(y,s )) = c(x\y)6{t - s ). (2) 
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Typical models have a quadratic form, g(x, y) oc (x - y) 2 for the decoherence term, corresponding 
to a fluctuating force which is independent of position. However, for a local bath one expects 
the correlation function to die off at some characteristic length scale (the correlation length of 
the environment), which has some important ramifications for decoherence. For a quadratic form 
of decoherence, the decay rate of the interference terms in the density matrix increases without 
bound, while for a local model the decay rate saturates at separations (between x and x') much 
larger than the correlation length of the environment, reflecting the independence of environment 
fluctuations at large separations. As it turns out, the quadratic form can be considered a short 
length scale approximation of a more detailed model. 

To consider the decoherence effects of an environment, simultaneous treatment of dissipation 
is necessary since decoherence and dissipative effects both generally arise from the same source 
(the interaction with a heat bath). For simplicity, I consider only linear dissipation, that is 

mx = p, p =-p + F. (3) 


As an example of quantum dissipative evolution, Dekker[7] has constructed a phenomenological 
master equation which includes ohmic dissipation and quadratic decoherence: 


dp 

dt 


il■+ 


(D X p T Dp X ) 
h 2 


[x, \p,p}] 


^rb> b> pW 


^r[x,[x,p}}. 


(4) 


The Caldeira-Leggett[8] master equation is obtained from a first-principle calculation for the effects 
of a simple thermal bath. With an appropriate choice of parameters for the Dekker model, the 
Calderia-Leggett master equation can be reproduced. 

Many open system models can reduce to the same classical phenomenology, particularly in the 
Markov regime, and yet have significant differences for a quantum sytem in that same regime. To 
illustrate this “richness” of quantum dissipative models, consider a rather generic oscillator bath 
model (following Zwanzig[9]): 


L = ]-mx 2 - U(x) + '52~yl < i> 


^liqp - <v(x)) 2 ]. 


(5) 


The classical calculations (the results of which are presumably reproduced in at least some limit of 
the quantum model) are relatively straightforward. The classical fluctuation-dissipation relation 
between the fluctuating forces and the nonlinear dissipation kernal emerges naturally, and in the 
usual Markov limit becomes: 

(fi{x,t)fj(y,s)) = k B Tr] tJ (x,y;t - s) =7] xj (x,y)2k B T6(t - s), (6) 


and a simple langevin equation can (at least in principle) be obtained: 

Xi(t) = ■ + /«(x(t), t) - 77 0 (x(i))xj. (7) 

For a homogeneous environment, the dissipation constant would be independent of position. 

Some observations about the Markov limit are in order. For the classical picture, the spatial 
correlations of the fluctuating forces are irrelevant. After all, the particle can only be in one place 
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at one time. For a quantum system one must consider superpositions between the particle at 
different locations, i.e. superpositions between different trajectories for the particle. My point is 
that different models may produce the same classical phenomenology, but have some important 
differences for the quantum case, in particular for the effective decoherence due to the environment. 

In order to help motivate some choices which will be required for the construction of the new 
model, consider a particle locally coupled to a scalar field. This particular model is a natural 
extension of one considered by Unruh and Zurek[10]. The action for this model is given by: 


L= 1 <Tr { i [*» - c’fV^) 2 ] + 6(t - x) - #(r, () - K(x)]} . 


( 8 ) 


This model produces approximately ohmic dissipation in one dimension[8,11]. In addition, one can 
extract from the influence functional the effective correlation function of the fluctuating forces[5, 
11 ]: 


< F ( x > r ) -F(y,s) >0= . f d d kk 2 ( — - 2 - cos (ut) cos(k • (x-y))l. (9) 

2{2n) J ( w ) 

This correlation function results from independent contributions from each mode of oscillation of 
the field. With some of the characteristics suggested by this local environment in mind, I now 
turn to the actual construction of the model. 

2 The Phenomenological Model 

The initial form of the evolution of the density operator is taken to be in the Lindblad(12] form 
(Schrodinger picture): 

% = L[p] = ±\H,p\ + il;[V>, V„'J + [ V„pV*\ = ±[H,p] + A L\p\, (10) 

for which there is a corresponding form for the Heisenberg picture V [O] which can readily be 
obtained from the cyclic properties of the trace. For a finite dimensional Hilbert space, this form 
is the most general for a completely positive dynamical semigroup. For infinite dimensional Hilbert 
spaces, it is a reasonable starting point. I will be focusing on the nonunitary part of the evolution, 
A L. 

The construction of the model is essentially the determination of the operators V^, subject to 
the constraint that the dissipation is ohmic (expressed as an operator condition). This constraint 
produces the “correct” classical phenomenology, but does not completely determine the model. 
However, linear dissipation almost forces the V h to be at most linear in momentum, that is 


V» = >M X ) - B m (x) • p. (11) 

Homogeneity and isotropy also serve to constrain the model. Assuming some sort of mode by 
mode interaction with a field, a reasonable choice is given by: 

(KJ = {a(k)e lk x — /?(A:)e ,k ' x k • p}. (12) 
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The discrete index p has been replaced by the continuous index k. The model is then completely 
specified by the complex functions a and 0. 

The resulting nonunitary contribution to the Schrodinger equation is given by the expression: 

AL[p] = - Jd t k i2 ^-(p - e’ kx pe“‘ kx ) - /d d A:^^(i{(k • p) 2 ; p} - e <kx k • ppk • pe' ik x ) 

_ . P) p}e“‘ kx ) - Jd d fc arn l Q ^/^ (e» kx [k • p,p]e“' kx ). (13) 


The position representation of the new model is given by: 

dp(x,x'-,t ) _ f f Jt _\ a ( k ) I 2 


dt 


Hamiltonian terms 


- (/-“ 


(1 - cos k(x - x*)) ) p(x , x'; t) 


- Q J dkRe[a (k)0(k)]k sin k{pc - x')j p(x, x'; t) 

- (i J dklm\a (k)0(k)]k sin k(x - x')j p(x,x';t) 

[jdk\m\w) **■*';*> 

{^J dk\0(k)\ 2 hk 2 cos k(x — x')^ ^ x > x '' t )• 


+ 

+ 


(14) 


The first nonhamiltonian term is responsible for decoherence. The corresponding noise spatial 
correlation is determined by a(k). The characteristic length should be on the order of the inverse 
of the “width” of |a(fc )| 2 in k space. The second nonhamiltonian term is responsible for the 
dissipation. Clearly the dissipation and other terms are more complicated in this new model. 
However, that would also be expected from a more elaborate first principle calculation. 

By examining the Eherenfest relations of physical quantities using L* , some interesting physical 
features of the new model emerge. By construction, the average position and momentum obey 
relations corresponding to ohmic dissipation: 


d{P) 

dt 


n m 


d{x) 

dt 


i {[H ' x]) 


(n 

m 


(15) 


where 


^ = fdk2Re(a (k)0(k))k 2 . 
m J 


(16) 


With only limited constraints on a and 0 (7 must be positive), the kinetic energy is seen to be 
thermalized: 

P?__ ksT, 

2m 2 


P 2 ksT 




(17) 


where 


7 


= 2^ - Id*l 

m J 


\m\\ u . k B T_l rjd Jc(k )I 


-hk\ 


d ’27 

The effective temperature is determined by a and 0. 


= - fd d k- 
7 J 


2m 


hk z . 


(18) 
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A low length scale approximation of the new model can be obtained by expanding the expo¬ 
nential terms in powers of k • x: 


A L[p] 




_ [k ■ x, |k • x, p]] - J<t d ^ l3t 2 h l k p ' l k ' p - '’ll 

|k x {k p p}|+ 


[k • x, [k • p,p]]. (19) 


The lowest non vanishing terms are second order, which exactly reproduces the Dekker master 
equation for 1 dimension. As a result, we can think of the Dekker or Caldeira-Leggett equations 
as a low length scale approximation for more general models. 

On the other hand, the Caldeira-Leggett master equation, 


A L[p] = ^[x, {p,p}\ - ^§i[x, [x,p}\, (20) 

can be considered a special case of the Dekker master equation, with the D xp terms equal to zero 
(which Dekker has argued should be the case) and an additional low momentum approximation 
which ignores the D xx term. With this type of special case in mind, we can construct a low 
momentum approximation for the new model which includes only the decoherence and dissipation 
terms: 


AL\p}= - fd d k&g£(p - e <k x pe- <k x ) - ({k • p, e <k x pe~ <k x }). (21) 

This would seem to be a likely starting point for applications of this model. However, this 
approximation is not a positive form for the dynamics. 

Finally, I would like to look at the Wigner representation of the new model, which has some 
interesting features. If we expand the terms of the evolution equation in powers of h (in the same 
manner as is typically done with the potential), 


W(q,p) 


1 + *L{y'( q )w) + V ^ 2 —y( 2 w+D d 2 "* 1 
mdq^ J ^dp [ Kq) * (2n)! am-i w 


+ A ±(pW) + 


+ Dr 


& 


I.'* 1 ” (/ 


n=l 


djpn+l 


w dp 2 

+ (O 


w + 


2n+l 


n=l 

^ +D ^m w+ ^ (hr+l {l 


J7 2R e(a(3)k 2n+1 \ d 2n+1 
(2n + 1)! ) dpW W 

[J dk {2n + 2)\J dp 2 "* 2 


+ D„^w + 52m 


2n— 


■(/ 


dk 


, Jm(a (3)k 2n+2 \ d <9 2n+1 
(2n + l)! ) dqd^ 

y ( h l d ± i w 

(2n)! )dp*\4d<l* V ) ’ 


( 22 ) 


the lowest order terms correspond exactly to the Wigner representation of the Dekker equation. 
The Wigner representation of the Dekker equation is a standard classical type of diffusion equa¬ 
tion. This illustrates the idea that the “classical” nature of the system emerges when coherent 
superpositions are not important in the dynamics. In this case, the relevant superpositions are 


423 



between different locations separated by distances on the order of the environment correlation 
length. 

The convolution theorem can also be used to write down the Wigner representation of the 
evolution: 


2 t j/Re[a* (£)#(£)] 

W(q,p) = (Hamiltonian terms) — — j dj/(p — j/)W ( q,p — j/)-- 


+ D ’ xi \h - rf w ( Q ’ p)+ J dp '~ 


\ 2 ,1 d 2 , (p-p , ) 2 WJ , / _ _ , 0 ^ 

~( t ^ + ±2 — W(q,p- p )• (23) 


Adq 2 


K 


One apparent effect in the new model is a spreading induced by these convolution terms. 

In summary, a new phenomenological master equation for ohmic dissipation and decoherence 
has been constructed which has completely positive dynamics. The new model has the features 
expected from local coupling to a homogeneous environment: specifically, the evolution is istropic 
and translationally invariant. Spatial correlations of the environment appear explicitly in the 
models. The new model also includes previous results as low length scale approximations. 


References 

[1] W. H. Zurek, S. Habib and J. P. Paz, Phys. Rev. Lett. 70, 1187 (1993). W. H. Zurek, Prog. 
Theor. Phys. 89, 281 (1993). 

[2] M Gell-Mann and J. B. Hartle, Phys. Rev. D 47, 3345 (1993). 

[3] M. R. Gallis and G. N. Fleming, Phys. Rev. A 42, 38 1990. 

[4] Some of these results have appeared in: M. R. Gallis, Phys. Rev. A 48, 1028 (1993). 

[5] M. R. Gallis, Phys. Rev. A 45, 47 (1992). 

[6] L. Diosi, Phys. Lett. A 112, 288 (1985). 

[7] H Dekker, Phys. Rep. 80, 1 (1981). 

[8] A. O. Caldeira and A. J. Leggett, Physica (Utrecht) 121 A, 587 (1983). 

[9] R. Zwanzig, J. Stat. Phys 9, 215 (1973). 

[10] W. G. Unruh and W. H. Zurek, Phys. Rev. D40, 1071 (1989). 

[11] M. R. Gallis, unpublished. 

[12] G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 


424 



